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Q_i, The process of proton emission from nuclei is studied by utilizing the two-potential 

approach of Gurvitz and Kalbermann in the context of the full many-body problem, 
r""". ; A time-dependent approach is used for calculating the decay width. Starting from 

an initial many-body quasi-stationary state, we employ the Feshbach projection op- 
erator approach and reduce the formalism to an effective one-body problem. We 
show that the decay width can be expressed in terms of a one-body matrix element 
If} \ multiplied by a normalization factor. We demonstrate that the traditional inter- 

pretation of this normalization as the square root of a spectroscopic factor is only 
valid for one particular choice of projection operator. This causes no problem for the 



o 



calculation of the decay width in a consistent microscopic approach, but it leads to 
ambiguities in the interpretation of experimental results. In particular, spectroscopic 
factors extracted from a comparison of the measured decay width with a calculated 
single-particle width may be affected. 

PACS numbers: 24.50. +g, 21.60.-n, 26.65.+t 

> 
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I. INTRODUCTION 

One of the classic problems in quantum mechanics is that of tunneling through a clas- 
sically forbidden region or, more specifically, the decay of a quasi-stationary state to the 
continuum. In nuclear physics, this manifests itself in the processes of a-decay in heavy 
nuclei and proton emission by proton drip-line nuclei. Of particular current interest are the 
lifetimes of proton emitters, especially in the lighter region of the nuclear chart, and the 
implications of this in nuclear astrophysics. 

Over the years, a number of different theoretical approaches have been used to describe 
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the decay process in nuclear physics, either by means of perturbation theory of decaying 
states or by time reverse study of resonance states via scattering theory 0, HI IE IH S El- 
Some authors solve the time-dependent problem while others use a stationary picture and 
make use of approximation methods such as the distorted-wave Born approximation or the 
semi-classical Wentzel-Kramers-Brillouin approach to evaluate the width Q. Other more 
accurate methods, such as R-matrix theory, are sometimes very sensitive to the channel 
radius giving dramatic variation in the calculated widths 0,0, 0,0]- The method of Gurvitz 
and Kalbermann 0, H,0|, also known as the two-potential approach (TPA) 0, 0] , is based 
on splitting the barrier potential into an interior and an exterior components. The inner 
potential binds the particle, which can then be described by a bound eigenstate of the 
relative Hamiltonian, while the outer potential acts as a perturbation that converts it into 
a quasi-stationary state (a wave packet), which can decay. 

An important shortcoming of all the above approaches, however, and in common with the 
descriptions of so many nuclear processes, is the approximate treatment of the many-body 
structure effects. In most descriptions of the proton-emission process the initial (74 + l)-body 
wave function is written as a product of an A-body wave function, describing the daughter 
nucleus, and the proton's single-particle wave function. The decay width is then written in 
the form of a single-particle width multiplied by a spectroscopic factor, which contains the 
many-body information of the system. This procedure, however, makes various assumptions 
about the relationship between the many-body problem and the effective one-body problem 
that have to be tested. In this work we consider the TPA of Gurvitz and Kalbermann and 
extend it to properly account for the many-body correlations. 

The standard reduction from a many-body problem to an effective one-body picture 
has been revisited in a recent study of radiative proton capture The work focused 

on one-body overlap functions and their associated equations of motion. The one-body 
overlap functions are obtained by integrating the product of the wave functions for an 
(A + l)-body system and its A-body subsystem over the coordinates of the latter. While 



the overlap functions are unambiguously defined, it was demonstrated in Ref. IllJ that 



useful 'auxiliary' one-body functions can be defined in several different ways. Natur ally the 



associated equations of motion differ for the three approaches considered in Ref. jllj. In 
the current work, we derive expressions for the proton decay width using two of the three 
approaches mentioned. The resulting decay widths have formally the same structure in both 
approaches, but the overall normalization factors differ. Only one of these normalization 
factors can be interpreted as the square root of a spectroscopic factor. This has consequences 
for the interpretation of experimental results and in particular for the determination of 
spectroscopic factors from decay widths. 

We start, in Sec. HH with the time-dependent Schrodinger equation and follow the stan- 
dard theory of decaying states [T3|. This method is briefly compared with the scattering 
approach to decay problems. In Sees. 11111 and II VI we use projection operator techniques and 
perturbation theory to derive an expression for the decay width in terms of the imaginary 
part of the pole in the Green's function matrix element. In Sec. 13 we describe the two al- 
ternative routes for reducing the many-body problem to an effective one-body problem. We 
compare the resulting two expressions for the width in Sec. IVII An alternative expression 
for the decay width is given in Sec. I VI II which shows more clearly the relation of the width 
to the spectroscopic factor. 
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II. PROTON EMISSION FORMALISM 



In this section we begin to recast the formalism of Gurvitz and Kalbermann in a form 
which is more convenient for our purposes. The approach of Refs. 0, E| starts with a square- 
integrable wave function, |^ ), which corresponds to the quasi-bound nucleus whose decay 
we are interested in. The initial wave function is close to the resonance state in the nuclear 
interior but decays rather than oscillates in the exterior region. This wave function cannot 
be an eigenstate of the full Hamiltonian or it would have a trivial time dependence and 
no decay would take place. Taking j^o) as the wave function at t = 0, we follow its time 
evolution using the time-dependent Schrodinger equation: 

ih^ t \m) = H\m)- a) 

This initial value problem is solved by using the one-sided Fourier transform (sometimes 
called a Laplace transform). We obtain 

-ih\1> ) + E$(E)) = H\i,{E)), (2) 

where $(E)) = f™dt e^ E+ie ^ h \i){t)) and e is a positive infinitesimal real number. Solving 
for \ip(E)) we have 

The probability amplitude for the nucleus remaining in the initial state after a time t is 
given by 



which is obtained by taking the matrix element of the previous equation with (i/j \ and 
carrying out the inverse Fourier transform. This is the Fourier transform of one particular 
matrix element of the many-body Green's function. In general the Lehmann representation 
of the latter contains contributions from many poles. However, the overlap (ipo\ip(t)) takes 
a simple form if the right-hand side of Eq. (@J) is dominated by the contribution of only 
one pole. In this case the decay rate can be extracted from the imaginary part of the pole 
location and we obtain a simple exponential function which describes the decay of the initial 
state. Thus the initial state \ipo) should be chosen to minimize the contributions from other 
poles. 

The same result can be found by considering the related scattering problem. A quasi- 
bound state can also be thought of as a resonance in the scattering amplitude. We consider 
scattering of a proton off the (A — l)-body system and we define a t-matrix, T, through the 
equation 

Ta\k)\^ A -i) = V\^ A ) } (5) 

where a\k) is the creation operator for a particle with momentum k = \J 2mE / h 2 and m 
is the reduced mass of the system. By standard techniques it can be shown that 

T = V + VGqT (6) 
= V + VGV, (7) 
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where G = l/(E — H + ie), G = 1/(E — H + ie), and H describes the free motion 
of the ejected particle with respect to the final state of the (A — l)-body system, i.e., 
H Q a^ (k)\ipA-i) = Ea){k) \1pA-1) and H = H + V. We observe that the poles of the t-matrix 
are given by the poles of V and G. Assuming V has no nearby pole we see that the poles of 
T are just those of the many-body Green's function. The width of the state is then given 
by the imaginary part of the pole location as in the previous case. 

In principle, the potential V is the sum of the interactions of the A-th particle with each of 
the particles in the (A — l)-body system. In nuclear physics this is generally approximated by 
a nucleon-nucleus optical potential. This approach, however, has the disadvantage of loosing 
track of the Pauli exchange correlations and other many-body effects. For bound states one 
of the principal effects is included through the use of the spectroscopic factor. However for 
scattering states there is, strictly speaking, no spectroscopic factor. Using Eq. (jU) as the 
starting point we can formally take the many-body effects into account while deriving an 
effective one-body equation. 



III. PROJECTION OPERATOR FORMALISM 



The main ingredient of Eq. (J3J) is the matrix element of the Green's function, M 



i^o] E _ H+ie \'4 J o)- Its expression can be simplified as done in Ref. [llj by using the projection 
operator formalism. We define a projection operator P = |Vo)(Vo| an d the complementary 
operator Q = 1 — P. The equation for the Green's function can be written as 

1 



(E-H) 



(E-H)G=\. 



(8) 



E-H + ie 

Acting on the left by P or Q and on the right by P we get the two equations 

P(E - H)(P + Q)GP = P, 

Q(E-H)(P + Q)GP = 0. (9) 
Solving the second equation for QGP and substituting this into the first equation we have 

1 



E - PHP - PHQ 



E - QHQ + ie 



QHP PGP = P. 



(10) 



Taking the matrix element of this equation with \ipo) and using the explicit form for P we 
find 



E- (Vol H + HQ 



E - QHQ + ie 



QH |-0o> 



(11) 



Since the left-hand side of this equation is a product of two terms, the poles of (^q\G\'^q) 
must coincide with the zeros of the multiplying factor. Thus the poles of the Green's function 
are given by 



E= (Vol ( H + HQ 



1 



E - QHQ + ie 
or, if we use a spectral representation of E _QXQ +ie , by 



QH |Vo) , 



E 



00 dE , \(MhQ\Ce>)\ 2 

E-E' + ie 



(12) 



(13) 
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where \(e') is the solution of the equation 

E\Ce) = QHQ\Ce), (14) 

normalized according to (Ce\Ce) = S(E' — E). Note that \iPq) is also an eigenstate of QHQ, 
but with energy E = 0. This state is excluded from the sum (integral) in Eq. ()13j) by the 
projection operator. In fact the only role of Q in this equation is to exclude the discrete 
state \ipo). The residue, R, of the pole of (i/)o\G\ipo) is given by 



R 



1 ~ lE^ HQ E-QHQ + U QH ^\ 



(15) 



This follows from a Taylor series expansion of Eq. (fTT|) in E. For the present problem 
R is close to one since we will choose |^o) to be an eigenvalue of H inside the nucleus 
and different from an eigenstate only in a region where the wave function is exponentially 
damped. Where |^o) is an eigenstate (ijjo\HQ is zero so the contribution from the second 
term will be exponentially small. 

More insight into Eq. ()12|) can be obtained by an alternative derivation. Consider the 
equation E = (i/;q\H\i/;e) / (i^oI^e) where \4>e) is an eigenstate of the full Hamiltonian, H. 
The wave function can be written as \i[)e) — {P + Q)\^e)- Following the standard Fesh- 
bach projection operator technique we can write Q\i/je) = Q E _QHQ +i£ QHP\ip E ). Using the 
explicit form of the projection operator P = |^o)(V'o| yields Eq. flT^j) immediately. 



IV. PERTURB ATIVE APPROXIMATION 

In general, Eq. ()12|) is highly non-linear and has many solutions. Every solution of this 
equation will give a pole of the Green's function, however, not every pole of the Green's 
function will necessarily be found by using this equation. For example, if has a defi- 
nite angular momentum only poles with that angular momentum can be found with that 
particular choice of |^ )- 

For the problem at hand, namely proton emission, we are not interested in the complete 
complexity of the Green's function. To obtain the decay width, only the imaginary part of 
the pole location of the nearby pole is relevant. We define an Hermitian Hamiltonian Hq 
such that H \ipo) = E \ip ) with H = H + 5H. Expanding the right hand side of Eq. (|T^|) 
about Eq and neglecting terms of order (E — Eq) 2 we obtain 

(M$H\*Po) + (HHQ^ n L, , ■ QH\ik) 

E Q - QHQ + it 



-(E — E ) 



^o\HQ——^——QH\^) 



(16) 



_dE" TU ' E - QHQ + ie Ir 'J e=Eq 

R ( (Vol^iVo) + (HHQ^ nnn^- Q H \^)] > ( 17 ) 

V Eq- QHQ + ie J 



where the last line has been obtained by solving for (E — Eq) in the first equation. This 
expression for the energy, Eq. (JT7J), contains a factor of R, the residue of pole of (iPq\G\^ ). 
Note that (iPq\HQ = (iPq\5HQ implies that R — 1 is of the order of 5H 2 . Since it multiplies 
a factor of order 5H, it will introduce terms of order higher than (E — Eq) 2 w 5H 2 and 
therefore it can be neglected. In the following, we thus take R — 1. 
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Let us emphasize that E is real since it is an the eigenvalue of the Hermitian operator 
H . The pole location, on the other hand, occurs at a complex energy. If the width of the 
state is large, for any reason not just proton emission, then E pa E does not hold (since 
E is real and the pole location has a large imaginary part) and the approximation fails. 
However, when the width is narrow then it will be possible to chose |Vo) or equivalently Hq 
such that perturbation theory is applicable. 

The equation for the energy is now linear in E and has the form typical of perturbation 
theory. The explicit connection with perturbation theory can be made when the unperturbed 
Hamiltonian is taken to be Hqq = Hq + Q5HQ, and not Hq as one might have expected. 
This definition is not only necessary in order to cast Eq. (JT3j) into the form of a perturbation 
expansion but at the same time eliminates the problem of non-compactness (H The 
non-compactness arises since Hq is chosen in such a manner that it does not go to zero at 
infinity. However both H and Hqq go to zero at infinity. 

The states \(e) occurring in the integral in Eq. (|T3Jl are eigenstates of H 00 , as is \ipo). In 
contrast to Eq. (|14j) the state |Vo) now occurs with energy Eq rather than 0. We can also 
rewrite (ip \HQ\^ E >) as (ifj \(H-H Q )Q\C E >) or (ip \(H - H 00 )Q\( E >) since (Vol is an eigenstate 
of both Hq and H 00 or, equivalently, since Q commutes with both Hq and Hqq. Using the 
previously defined 5H = H — H Q , Eq. (|T3j) can be rewritten as 

g ^ +(tfgw+ r tf »)f (18) 

J-oo E - E> + ie 

^ 0+(Wo>+ r^ ^jjM , (19) 

J-oo Eq - E +l€ 

where 



5H' = H- H 00 = P5H + 5HP - P5HP. (20) 

Equation (|18|) for E is not in the form of perturbation theory since \( E ) is an eigenfunction of 
H Q0 and not Hq. However, Eq. (J19)) is. The only approximation made in deriving Eqs. (|18|) 
and (fTTJj) was the replacement of E in the denominator on the right hand side by Eq, i.e., 
we assumed that perturbation theory is valid. 

We emphasize that neither | Vo) nor \Ce) is an eigenstate of the full Hamiltonian H and 
our derivation depends crucially on this point. Instead, |Vo) is an eigenstate of Hq and 
Hqq while \( E ) is an eigenstate of QHQ and Hqq. The fact that both wave functions are 
eigenstates of H 00 allows us to interpret the right-hand side of Eq. (J 19)) as the first few terms 
of a perturbative expansion. 



V. CHOICE OF H, 



The proper choice of H was discussed in a very convincing and direct way by Gurvitz 
and Kalbermann 0, for the problem of a single particle in a local potential well. The 
reader is referred to those papers for specific details. The basic idea is to take an H such 
that at infinity the potential goes to a finite value, larger than E , rather than to zero. This 
causes the state under consideration to be bound, but if the decay width is small this new 
state should be very close to the scattering state. We now generalize this concept to the 
case of A interacting particles. 
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We consider the case where there is one open channel. Asymptotically the wave function 
describes a free proton and the bound (A — l)-body system. Thus we take 

6H = -J dr a\r)\^ A ^)V(r)(^ A ^\a^), (21) 

where |$a-i) is the ground state of the (A — l)-body system, r is the relative coordinate of 
the proton and the (A — l)-body system. The integral covers the whole space. If V(r) is 
taken to be greater than Eq outside the range of the nuclear potential, ro, it will prevent the 
initial state \ipo) from decaying. At the same time, we want the perturbing potential V(r) 
to be zero inside the nucleus so that it does not modify the wave function in the interior. 
With this ansatz expression (|T5j) for E becomes 

«(,)„(,)*,<,) + r M IJ*«WyWI 3 , (22) 

J J-oc E - E' + ie 

where 0o( r ) — (^A-i\0'(' r )\4'o) is a spectroscopic amplitude (i.e. a one-body overlap function 
involving bound many-body states) and 4>fAt) = ($A-i\ a ( r )\(E') is an optical model wave 
function (cf. the discussion in Refs. flflL llll])- The projection operator Q does not have to 
be included explicitly since the integral does not include the discrete state 4>o(r). 

A miracle has occurred here. Due to the choice of 6H, which is physically motivated, 
the expression has reduced to an effective one-body problem. All the many-body aspects of 
the problem are contained in the one-body overlap functions, <po(r) and <pE'( r )- We stress 
that the only approximation made so far is that second order perturbation theory was used 
to justify replacing E with Eq in Eq. JTJJ). For states with narrow widths this should be 
acceptable. We have not assumed that \ipo) is a product state or made any other assumptions 
regarding its structure. 

Equation (|T3j) for \( E ) can be written as E\( E ) = H 00 \( E ). In the neighborhood of E — E 
the Hamiltonian in this equation for \Ce') can be approximated as 

SH' « - J dva\v)\^ A ^)V m {v)(^ A ^\a{v). (23) 

For an appropriately chosen Voo(r), the approximation made here is the same as that of 
Eq. (2.15) in Ref. The basic argument given there proceeds as follows: the eigenfunction 
of H at resonance will be large when the AtYt particle is in the nuclear interior. However the 
state \(e>) will be small due to the projection operators in QHQ. It will look like the real 
state in exterior region but be suppressed in the interior. The form of Voo(r) will ensure this 
if we take it to be in the exterior region and repulsive (greater than Eq) in the interior 
region (see Fig. 3 in Ref. [H). This is the opposite of what we did for V(r) which was large 
in the exterior and in the interior. 

Despite its nice form, Eq. (|2*2*|) is not useful until we specify how to calculate the functions 
4>q (r) and 0£/(r). Following Ref. [ll| we set up one-body equations for these functions using 
the Feshbach projection operator technique. We start with the exact scattering state 1$^) 
and write the equation-of-motion for the corresponding one-body overlap function (which in 
this case is the Feshbach optical model wave function) and see how it is modified when the 
perturbing Hamiltonian SH is added. The one-body overlap function for the exact scattering 
state is given by (Tl|, Q 03 



£0(r) = j dr'dv"(<$> A „ 1 \a(v) \H + HQ F E _ ^ HQ J Q F Il\ a) {r")\<$> a „ x ) 
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Af (r", rT V(r') (24) 
dr'H(r,r')4>(r'), (25) 

where Pp = J drdr' a^(r)|$ J 4_ 1 )A/'(r, r')~ 1 (^A-i\ a ( r '), Qf = 1— -Pf and the norm operator is 
jV(r,r') = ($A-i|a(r)a t (r / )|$ A _i). Since <5#Qf = Qf^# = 0, it follows that <p (r) satisfies 
the equation 

E<f) (r) = J dr' [H(r, r') + Af(r, r')V(r')\ o (r'). (26) 

Note the explicit asymmetry of this Hamiltonian. Even in the absence of inelastic processes 
it is not Hermitian. An analogous potential in alpha emission has been strongly criticized 
[TH Il7l ] in the st udy of alpha emission. Two methods of addressing the problem were 
discussed in Ref. |ljj. The first is to define a new amplitude 0(r) = J dr' A/"(r, r')~ 1 / 2 0(r / ) 
and take 

5H= -J drdr'drV^^A/^y^^ (27) 

The equation for 0(r) is then 

E0 o (r) = J dr' [H(r, r') + 8(r, r')V(r')} o (r'), (28) 

where 7Y(r, r') is given in Ref. 0] as 

H(r, r') = J dr"dr"'M(r, r")-^ 2 n(r", r"')Af(r"', r') 1/2 . (29) 

A similar development holds for </>e(i") and we obtain the corresponding equation 

EMr) = J dr' [H(r, r') + 6(r, r')V 00 (r')} Mr'). (30) 

The equations we find for 0o( r ) and 4>e{^) are quite remarkable. As previously noted, we 
have reduced the many-body problem to an effective one-body equation in which all many- 
body effects are contained in 7i(r, r'). The perturbing potential is taken to be local in both 
cases. This is not really an additional approximation since the justifications given previously 
still hold. 

Another way of getting a symmetric form for the additional potential is by using a 
projection operator given in Ref. 



P F = J dr [o(r) +at(r)]|$ j4 _ 1 )($ A _ 1 |[o(r) +o+(r)] , (31) 

for which Eq. (|2*Bj) reduces to the Dyson equation of many-body theory. For time reversal 
invariant states the associated norm operator reduces to the unity operator. Repeating the 
derivation with this new projection operator we obtain 

£0o(r) = f dv'[H M (r,r') + 5(r, r')V(r')]cf> ^l (32) 

E(j> E {r) = J dr'[n M (r,r') + 5(r, r')V 00 (r')] <j> E (f), (33) 
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where ?^(r, r') is the mass operator that occurs in the particle-hole Green's function 
In general there is no simple relation between TCj^(r, r') and 7Y(r, r'). The equations for 0(r) 
and 0(r) are now formally the same. At this point we would expect the widths calculated 
with these functions to be numerically similar. The only difference in the two treatments is 
the form of V(r) and Voo(r). When the argument leading to the justification of perturbation 
theory is correct we expect both approaches to work equally well. 

The argument for Eq. (j23j) can be restated using the one-body functions defined above. 
The resonant state, that is the solution of the exact Schrodinger-like equation Eo<f)Q Ue (r) = 
J dr' 7i(r, r')(pQ nc (r'), grows exponentially as we go into the classical forbidden region under 
the Coulomb barrier from the outside and the wave function is exponentially large in the 
interior. On the other hand, <fig (r) coming from \Ce), the solution of the many body-state 
defined with the Hamiltonian QHQ, is exponentially suppressed due to the orthogonality of 
l^o) an d \(e)- If the suppression is large (i.e. the state under consideration is narrow) the 
precise form of \(e) in the interior is not important as it is essentially zero. The approximate 
form of H 00 is chosen so the approximate \(e) is also zero in the interior. As we will see 
later we do not need in the interior but only outside at some radius, r . 



VI. TWO EXPRESSIONS FOR THE WIDTH 



The decay width is given in terms of the imaginary part of the pole location in Eq. 
We find 

2 



or 



T « 2tt 



2tt 



dr <j>* (r)V(r)<f> Eo (r) 



dr 0*(r)V(r)0 Bo (r) 



(34) 



(35) 



depending on whether we use the bar-ed amplitudes or the standard overlap functions. As 
discussed below, we expect the approximations leading to the two equations to be valid 
simultaneously, so that the resulting two expressions for the width will agree. We stress 
again that the reduction to an effective one-body problem is not an approximation but 
emerges rather naturally from the formalism. In a pure one-body problem 0o( r ) and 0o( r ) 
would be normalized to unity whereas here they are normalized to the factors v^o and \/So, 
respectively. The normalization factors contain the many-body aspects of the problem. For 
the standard overlap function <po(r), this normalization is the well-known spectroscopic factor 

(see also Ref. [HI)- It is useful to define normalized one-body functions 0o( r ) an d 



r via 



<p {r) = a 
functions 



/ 3o0o(r) and 
Eo (r) and q 



E { 



i.e. 



jdr 



fdi 



The wave- 



rj = V^o<Pol 

r) describe scattering states and are normalized asymptotically 
at large r. If we assume the one-body Hamiltonians TC and ft to be local (more will be 
said about this later) we can take over verbatim the development given in Refs. 0, H, 0| . 
The only difference between the treatment presented there and the approach shown here is 
the presence of the factors S (for the standard overlap functions) and S (for the bar-ed 
function). To make the dependence on these normalization factors explicit we rewrite the 
last equations as 

2 



2nS n 



dr 0S(r)K(r)0 Eo (r) 



(36) 
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and 



2nS 



(37) 



To further understand the results we peruse one particular line of development. As 
previously noted, we take V(r) to be zero inside the nucleus, V(r) = for r < r . For 
t > t , we have J dr'W^(r, r'VA^(r') = Ec/)e(t) since Voo(r) vanishes for large r. For 
definiteness we follow Refs. P, |9J] and take tq to the maximum of the potential. For the 
example of fluorine discussed below, this corresponds to about 5 fm. Thus we write the 
integral in Eq. as 



dr9{r-T Q )4>* Q {v)V{v)4> EQ {v) (38) 
drdr' 9(r - r )$(r) {H M (r,r') + [-E + V(r)]5(r, r')} ^ (r), (39) 



with a similar equation holding for the bar-ed quantity. Naively one might expect this to 
vanish and indeed it would if the step function 9 were not there. In the following we shall 
consider the case in which the non- locality in TCj^(r, r') is small enough to be properly 
described by a local potential and an effective mass m^r). In this case one has 



n M (r,r>) 



2mfc(r 



5(r,r') 



(40) 



In this equation Hm is correctly evaluated at the energy E of the scattered state, therefore 
mfc(r) accounts for the sole spatial non-locality and differs from the usual definition of 
effective mass m*(r) ET 



20j ] . When Eq. (|40j) holds, the integral, Eq. (J39)) . can be simplified 
through integration by parts and, for spherically symmetric potentials V(r), Vm(t) an d 
effective mass rrik(r), one obtains 



/ = 



h 2 

[0r*o( r o)0ri?o( r o) - 0r£o( r o)0r*o( r o)] 



2/7lfc 

h 2 



2fhk ( ro ) ^r*o( r o)0ri?o( r o) - <PrE {r )<j)'*o( r o)] , 



(41) 

(42) 



where the wave functions are assumed to factorize as 0(r) = ^ £ ^-Y[ n (6, ip) and Yf l (8, if) are 
spherical harmonic functions of the angular variables 9, (p. The prime denotes the derivation 
with respect to the radial coordinate r. 

In deriving these last two equations we made an additional approximation, namely that 
the potentials are local in the vicinity of ro [or at least the non-locality is restricted to an effec- 
tive mass mfc(r)]. The two equations will be simultaneously valid only if the norm operator is 
local, i.e. J\f(r, r') = 7V(r)<5(r — r'), in the vicinity of r . In this case fhk{r ) / rrik(r ) =Af(r ) 
and the two equations are identical. Since r is in the tail of the density distribution, one 
expects the norm operator to be unity in its vicinity, which means that the two expressions 
for the width should give the same result for realistic models. This is verified numerically in 
Fig. Q where the function <p T (r) (solid lines) and 4> T (r) (dashed lines) are displayed for the 
3/2~ and 3/2 + decaying states of 17 F at respective excitation energies of 4.64 and 5.00 MeV 
(the proton threshold is at 0.60 MeV). The functions are calculated with the self-consistent 
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Green's function method of Refs. 21. HHJ. The 3/2+ state is well explained by the nuclear 
mean field approach and provides a typical example of a strong state, for which Sq ~ 5n. 
The 3/2~ state, on the other hand, is a typical example of a weak state for which Ref. [23j 
gives So ~ 0.04 and S ~ 0.14. It is worth to emphasize the difference between these two 
cases. For strong states the quenching of the spectroscopic factor is due both to short-range 
correlations and to the coupling to other excitations of the system {22|. Nevertheless they 
maintain a strong single-particle structure and the orbital occupancy is of the order of unity. 
Weak states, instead, have a more complicated structure and can be seen as collective ex- 
citations on their own rather than having a single-particle character. As a consequence the 
one-body spectroscopic factor can be one order of magnitude smaller or less. In this case 
the functions 0(r) and 0(r) sample the one-body substructure in a different way, whence 
the difference in their normalizations with S larger than S [11]. For 16 + p, the nuclear 
interaction becomes negligible beyond 5 fm typically and Fig. ^ shows that 0(r) and 0(r) 
are equal beyond this radius, which means that the norm operator is unity. In principle, 
there is no problem calculating the width from a microscopic model: one may define two 
different functions 0(r) and (f)(r), which have two different normalization factors a/S 7 and 
\/S, respectively, but the estimate for the width will be the same with both of them. In fact, 
a similar formula would be valid for any amplitude of the form <f) u (r) = J dr' 7V(r, r') !/ 0(r') 
for arbitrary since jV(r, r') is unity around tq. 

Let us emphasize that for the actual calculation of the width by the above formulas, a 
microscopic model based on the harmonic oscillator basis, as in Refs. 0, 23[ may not be 



the best choice since the wave function is not expected to be very precise in the vicinity 
of r . This is demonstrated in Fig. ^ where the 3/2 + resonant wave function of the phe- 
nomenological mean field Woods-Saxon potential of Ref. [24] is shown for comparison. The 
agreement with the microscopic wave function is good in the interior but deteriorates above 
4 fm (where the approximation of the potential by a harmonic oscillator breaks down, see 
Fig. |2J). The difference between both models is particularly large in this case since the 3/2 + 
state is wide (Fo = 1.5 MeV); for narrow states, the harmonic oscillator approximation 
could be sufficient, a possibility which will be explored elsewhere. Let us finally stress that 
the phenomenological mean field potential of Ref. |24[ does not reproduce the 3/2— state 
because its structure is not well approximated by an 16 core plus a proton. 

Equations (J36|) and (J37j) indicate that there may be a serious problem in extracting 
spectroscopic factors from measured decay widths. The standard method for determining 
a spectroscopic factor involves dividing an experimental width by a single-particle width 
calculated with a phenomenological model. However, it is not clear a priori whether the 

phenomenological wave functions are good approximations to 0o( r ) and 4>e{y), or to o ( r ) 
and 4>e{ t )', hence, it is not obvious whether dividing the experimental width by the result 
of a single-particle calculation provides the spectroscopic factor S or the normalization S 
of the auxiliary function 0(r) (or the norm of yet another one-body function). Normally, 
one assumes in proton emission studies that 0o( r ) and </>_e(i") can be equated with the wave 
functions obtained from phenomenological potentials (see for example Ref. 0]). On the 
other hand, in the context of some alpha emission studies it has been argued very strongly 
that 4> (r) and 4>e (r) correspond to phenomenological wave functions 3, ^} . If the latter 



is true then the experiments would be sensitive to Sq rather than to Sq. For strong states, 
with a clear core-plus-particle structure, this is mainly a philosophical issue, since S ~ S 
holds. For weak states, however, Sq and S can be significantly different from each other. 
We have attempted to resolve the ambiguity outlined above by calculating the effective 
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r [fm] 



FIG. 1: Radial part of the one-body overlap functions 4> (solid lines) and of the auxiliary functions 
4> (dashed lines) for the lowest 3/2~ and 3/2 + states of 17 F, calculated via the self-consistent Green's 
function method. For the 3/2 + state, both functions are nearly identical, and a phenomenological 
wave function, with a normalization chosen to be close to that of the microscopic wave functions, 
is shown for comparison (dash-dotted line). 



local potentials corresponding to </>(r) and <j)(r) and comparing them with typical phenomeno- 
logical potentials. This is done by inversion of the local (radial) Schrodinger equation: 

where the effective potential V e g(r) is the sum of the interaction potential (nuclear + 
Coulomb) and the centrifugal term. The effective potentials corresponding to the radial 
wave functions of Fig. ^ are shown in Fig. El For the strong 3/2 + state [1 = 1 centrifu- 
gal term) the three potentials are in reasonable agreement, except above 4 fm where the 
potential extracted from the microscopic functions asymptotically approaches the harmonic 
oscillator potential that generated them. For the weak 3/2~ state [I = 1 centrifugal term) 
the potentials deduced from 0(r) and 0(r) display a singularity and are very different from 
traditional phenomenological potentials. This singularity occurs because the zeroes of r (r) 
and 0" ( r ) occur at different radii. This is probably not an artifact of the model but a real 
effect and arises from the relative sign of the and 2ftu contributions to the spectroscopic 
amplitude. Let us remark that Eq. (|43J) assumes a constant effective mass rrik(r) = m. We 
have checked that introducing a realistic effective mass does not lead to a simpler potential 
for the weak state; more will be said about this elsewhere. Let us finally emphasize that, 
though their energies are close to one another, the strong and weak states correspond to very 
different potentials, which suggests that a strong energy dependence is a necessary feature 
for such potentials (see the discussion in the next section). In conclusion, this example shows 
that constructing reliable phenomenological local potentials for extracting spectroscopic fac- 
tors from experimental cross sections is non-trivial. Moreover, since the characteristics of 
the potential depend strongly on the state it is difficult to determine which normalization 
(5o or <So) would be extracted from a comparison with the experimental data. 
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FIG. 2: Effective (interaction + centrifugal) potentials corresponding to the wave functions of 
Fig. ^ The 3/2~ potentials have an I = 1 centrifugal term and a singularity around 3-4 fm; the 
3/2 + potentials have an I = 2 centrifugal term and are regular. The dash-dotted line indicates the 
potential corresponding to the phenomenological 3/2 + wave function. 
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VII. ALTERNATIVE EXPRESSIONS FOR THE WIDTH 

Let us now return to the expression for the decay width and establish a link with known 
results. In Refs. p, la l| the spatial derivatives were evaluated using a special form for V(r) 
and Voo(V)- m those references these auxiliary potentials were chosen such that the total 
potential for rO (r]was constant outside a given radius and the one for r £ o (r) was constant 
inside. In Refs. [| || the separation radius is at the maximum of the potential barrier, 
whereas in Ref. [5] the advantage of using a radius of the order of the nuclear range (our r ) 
is pointed out. For this last case, the width is then given by 



fc2 "I 2 

Trot 



2m k (r ) 



|0ro(r o )<Wr o )| 2 , (44) 



where a = \j2m k (r )[V(r ) - E ]/h 2 . 

An alternative approach is to exploit the Wronskian form of Eqs. (j41|) or ()42j) . If V{r) 
is zero for radii less than the outer turning point and Voo(r) zero for radii greater than 
the inner turning point then there is a region where r o( r ) and r _B o ( r ) satisfy the same 
differential equation. If the potential is local the Wronskian is a constant. Assuming a thick 
barrier, there will be a region where r o( r ) is an irregular Coulomb function, the regular 
Coulomb function having decayed away, and where 4> T E ( r ) is a regular Coulomb function, 
the irregular Coulomb function having decayed away. All that is required is to determine 
the proportionality constants. These have a simple expression for a constant effective mass 
m k (r) = m: for r £ o ( r ) we have r _B o ( r ) = ^/2m/h 2 k TrF(k r), with k = ^/2m£ /ft 2 , while 
for (f) r o( r ) the proportionality constant can be written as 4> T o( r o) /G(k r ). Here F(k r) and 
G(k r) are the regular and irregular Coulomb functions respectively, the Wronskian of which 
equals — k . The resulting width is then, as in Ref. 



!rkn \Mr )/G(k Q r )\ 2 . (45) 



m 



This can be simplified further. The wave function <p r o(r) is normalized to the appropriate 
spectroscopic factor, S^. We can also consider the true scattering wave function at resonance, 
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l true I 



In the interior it will behave like 4> r o( r ) while in the exterior region it will behave 



like G(k r). Normalizing it to G(k r) in the exterior region we obtain for the width: 



r r 



S ftv 



£ dr\<f>^(r)\ 



(46) 



where vq = hko/m is the asymptotic velocity. The exact value of the upper limit on the 
integral is not crucial and we take it to be the outer turning point (see numerical justification 
below). 

Equation (j4*o^l can also be derived in a more transparent way (see also Refs. [lH litSI]). 
Consider the Schrodinger-like equation \H(E) — E]<p(r,E) = for the overlap function. 
Since the effective one-body equation we have been considering can depend on the energy 
we keep an explicit energy dependence in the Hamiltonian. Differentiating this equation 
with respect to E we find 



[H(E) - E] 



d<p(r,E) 
dE 



1 - 



dH(E) 
dE 



(r,E). 



(47) 



Next we multiply by <f>*(r, E) and integrate up to some radius r\. If H(E) is not hermitian 
then 4>*(r, E) should be replaced by the time reversed state. If the potential is local in the 
vicinity of r; we can integrate by parts on the left hand side. Assuming spherical symmetry, 
this gives us 



h 2 



2m k (ri) 



dE 



dr(f)*{r,E) 



dE 
dH v (E) 
dE 



(r,E), 



(48) 



where the prime denotes a partial derivative with respect to r and H T (E) is the radial 
Hamiltonian. If we take 77 to be outside the range of nuclear force, (f) z (r, E) can be written 



as 



(r,E) « cos 5 (E) F(kr) + sin 5 (E)G(kr). 

r>ro 



(49) 



At a narrow resonance the phase shift 5(E) will be rapidly varying so we expect that the 
largest part of the energy dependence will come from the phase shift and not from the 
Coulomb functions. The energy dependence of the Coulomb functions will be minimized 
if the radius is chosen to be near the outside turning point. For example, for large r the 
regular Coulomb function will have a sin(kr) dependence. Differentiating with respect to E 
will give r(dk/dE) cos(kr), which diverges for large r. As r decreases to the turning point 
this asymptotic form for the wave function breaks down. However a similar argument using 
exponentials holds inside the turning points. Thus near the outside turning point we have 



h t (r,E) _ d5(E) 
~dE ~ dE 



- sm 5(E)F(kr) + cos 5(E)G(kr)}. 



(50) 



We have checked this relation numerically for Woods-Saxon plus Coulomb potentials and 
verified that for widths less then 15 KeV the error does not exceed 3%. That the radius 
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should be chosen near the outer turning point was also confirmed numerically. The Wron- 
skian relation can now be written as 



h 2 k d6(E) 
2m k (r t ) dE 



dr(p*(r, E) 



dH r (E) 
dE 



(r,E), 



(51) 



where we have taken ri = r t to be the outside turning point radius. 

At the resonant energy we expect the energy variation of the phase shift to be a maximum 

so the resonance energy occurs when J Q n drcf)*(r, E) 1 — 9H q^ r ( r ; E) is a maximum. At 

the resonance energy d5(E )/dE = 2/T , and we have 



/ o rt dr0*(r,£ o ) 



dH r (E ) 
dE 



(r,E ) 



(52) 



where now v = hk /m is the asymptotic value of the velocity. For bound states the Spee- 



ding) 
dE 



the overlap func- 



troscopic factor can be written as So = [J °° dr</>*(r)</> r (r)] / | J °° dr<j)*(r) 

(see Ref. (I?!), which does not depend on the specific normalization o 
tion 0(r). By extending this relation to define the spectroscopic factor for resonant states, 
we recover Eq. flJBJ). 

We finally note that the expression (J52j) of the width is independent of the choice of 
r (r, E) or r (r, E). Its denominator can be rewritten as 



D 



dr<f)*(r,E) 



_d_ 

dE 



[E — H X {E)\ \ <f) T (r, E) 



(53) 



The bar-ed amplitudes are defined as 0(r, E) = J dr A/"(r, r') _1//2 0(r, E) and a corresponding 
expression for the Hamiltonian is 8(E) = A/" 1 / 2 \H(E) — E] M 1 ^ 2 +E. Thus the denominator 
is invariant under this transformation, as well as under a general transformation with an 
arbitrary power of A/"(r, r'). All that matters is that H T (E) and 4> T (r, E) are consistent with 
one another. These wave functions are phase equivalent and any of them can be used. 



VIII. DISCUSSION 

We have embedded the elegant Gurvitz-Kalbermann approach of proton emission 0, H, 0| 
into a full many-body picture. We have reduced the formalism to an effective one-body prob- 
lem and demonstrated that the decay width can be expressed in terms of a one-body matrix 
element multiplied by a normalization factor. At first sight, this result agrees with the 
standard procedure for extracting spectroscopic factors from measurements via dividing an 
experimental width by a calculated single-particle width (see, for example, Ref. |7(). The 
present work, however, clearly demonstrates that this procedure for determining spectro- 
scopic factors is only valid if the phenomenological potential used to generate the single- 
particle width corresponds to the potential in TCm [ see Eqs. (|3~2I) and (|3*3*|) ]. It is not a priori 
clear that this is actually the case. In fact, the authors of Refs. |l6l.ll7| (and prior to that the 
authors of Ref. [28| ) have argued strongly that the phenomenological potential approximates 
the potential in H [see Eq. While the studies of Refs. (ly, |l7| were carried out for 

alpha decay, the arguments given there can be carried over to a description of the proton 
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emission process. Furthermore, Eq. suggests that JP dr<fi*(r) 1 — 9H q^ <j>r(r) is the 

appropriate observable that can be extracted from proton emission experiments. 

Besides the ambiguity regarding whether the standard spectroscopic factor or an auxiliary 
normalization is extracted from the experimental procedure, it has been demonstrated that 
constructing a reliable phenomenological potential is non-trivial. The situation is quite 
complicated since the interaction with the nuclear medium strongly depends on the initial 
state of the ejected proton. For states with a clear core-plus-particle structure (i.e., with 
a large spectroscopic factor), traditional phenomenological potentials seem to provide good 
approximations to the nuclear mean field and spectroscopic factors can be determined from 
proton emission studies. In this case, the spectroscopic factor extracted from the experiment 
can be safely compared to the results of nuclear many-body calculations (note also that for 
large So values, So and So are approximately equal [ll| and the distinction between the two 
approaches discussed here becomes irrelevant). 

For weak states, which have a more complicated many-body structure, standard phe- 
nomenological potentials do not give a proper approximation to the nuclear medium, as 
shown by the radial shape of the 3/2~ states in Fig. |21 Also, as discussed in Ref. fill ], the 
dependence of the spectroscopic factor on the energy derivative of the effective one-body 
Hamiltonian implies that the nuclear medium must be strongly energy dependent. This 
feature, which is missing in most phenomenological optical potentials, is also confirmed by 
the numerical results displayed in Fig. Thus for weak states simple potential models are 
probably not valid for either 0(r) or 4>(r). States that are neither weak nor strong will also 
have to be dealt with on a case by case basis. 

Since 0(r) and 0(r) are identical for large radii they are phase equivalent and elastic 
scattering experiments cannot distinguish between them. We conclude that additional ex- 
perimental input, together with an accurate derivation of the optical potential based on first 
principles, is required in order to resolve the question regarding which one-body Hamiltonian 
is most appropriately approximated by a phenomenological model. 
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